The phase diagram of a gauge theory with fermionic baryons 
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The fermion-sign problem at finite density is a persisting challenge for Monte-Carlo simulations. 
Theories that do not have a sign problem can provide valuable guidance and insight for physically 
more relevant ones that do. Replacing the gauge group SU(3) of QCD by the exceptional group G2, 
for example, leads to such a theory. It has mesons as well as bosonic and fermionic baryons, and 
shares many features with QCD. This makes the G2 gauge theory ideally suited to study general 
properties of dense, strongly-interacting matter, including baryonic and nuclear Fermi pressure 
effects. Here we present the first-ever results from lattice simulations of G2 QCD with dynamical 
fermions, providing a first explorative look at the phase diagram of this QCD-like theory at finite 
temperature and baryon chemical potential. 

PACS numbers: ll.30.Rd 12.38.Aw 12.38.Gc 12.38.Mh 21.65. Qr 



Finite fermion density continues to be a serious chal- 
lenge for Monte-Carlo simulations due to the fermion- 
sign problem [l|, y] . The sign problem appears in many 
areas of physics, but is of notorious importance to dense 
quark systems, especially in nuclei, heavy-ion collisions, 
and compact stellar objects. An alternative are models 
and continuum methods which do not have this type of 
problem [3|-l7|. However, these usually require approx- 
imations, and cross checks through lattice simulations 
remain desirable to improve systematic reliability. 

To provide support from numerical simulations, two 
major strategies have been followed. One is to replace 
the baryon chemical potential by some quan tity more 
amenable to simulations, e.g. imaginary [8|-[l0| or isospin 
prl . [l3 | chemical potential. The other is to replace the 
theory with one accessible through numerical simulations 
at finite density. However, such theories usually differ 
from the original one in more or less important aspects. 

One very well studied replacement of QCD for strongly 
interacting matter at finite density is two-color QCD |13l - 
[l7| . In this case, the baryons are bosons instead of 
fermions, however. This leads to profound differences, 
such as Bose-Einstein condensation of a baryon super- 
fluid with a BEC-BCS crossover at high densities in- 
stead of the usual liquid-gas transition of nuclear matter. 
While two-color QCD has many interesting aspects that 
deserve to be studied in their own right, the quantum ef- 
fects due to the fermionic nature of baryons are expected 
to play a very significant role for nuclear matter and es- 
pecially in the physics of compact stellar objects jlSf . 

Therefore, a more realistic replacement theory in this 
regard should contain fermionic baryons. One possibility 
is the strong-coupling limit [l9[. In order to maintain 
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the connection with the continuum, however, we employ 
here a different theory without sign problem for Monte- 
Carlo simulations. It is obtained by replacing the SU(3) 
gauge group of QCD with the gauge group G2 [20|- All 
color representations of this theory are equivalent to real 
ones. As a consequence the Dirac operator has an anti- 
unitary symplectic symmetry which for Nf flavors leads 
to an extended Pauli-Giirsey SU(2A^/) x Z2 flavor symme- 
try [13|,|20| , and is thus expected to have a non-anomalous 
component even for a single flavor. In this paper we study 
the phase diagram for a single Dirac flavor of Wilson 
fermions, corresponding to a continuum SU(2)xZ2 ex- 
tended flavor symmetry in the chiral limit. Spontaneous 
or explicit breaking reduces this to U(l)xZ2 [2ll]. The 
unbroken U(l) relates to the baryon number to which the 
baryon chemical potential is coupled. The anti-unitary 
symplectic symmetry of the Dirac operator in a real color 
representation implies a two-fold degeneracy of its eigen- 
values and hence that the fermion determinant remains 
positive at finite baryon chemical potential even for a 
single flavor [ll,[il. 

The physical bound states of this theory, besides the 
usual quark-antiquark and three-quark states, also con- 
tain hybrids of one quark with three gluons, as well as 
diquarks and further bound states with more than three 
quarks [201 . Thus, the hadronic spectrum contains both 
fermionic as well as bosonic baryons and mesons. These 
bound states are created by very similar interactions as 
in QCD, i. e. by a potential which rises linearly with 
the separation of the quarks before string-breaking sets 
in [20I. l22 - l25l |. In fact, in a simple quark model picture 
with light quarks, one has a baryon-mass hierarchy where 
for a sufficiently small quark mass nucleonic three-quark 
states will be the second lightest baryons, above the light 
diquarks as the would-be- Goldstone bosons of the ex- 
tended flavor-symmetry breaking. Hybrids, tetraquarks 
and other baryonic multi-quark states are expected to be 
much heavier. 

Another appealing aspect of the G2 gauge group is 
that it might be deformable to ordinary QCD by break- 
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FIG. 1: The left panel shows an example for the mass correla- 
tor fits performed. Since in most cases no clear mass plateaus 
appear due to closeness of the excited states, conservative er- 
ror estimates have been used, which are shown in the example 
as error bands. The right panel shows the masses obtained in 
this way. Note the shifted /3 values to disentangle the plot. 
All results are for k = 0.15625. Only the range j3 = 0.9 to 
P = 1.0 has been used in the phase diagram calculations. 

ing G2 down to SU(3) via a Higgs mechanism [20, [26| . 
although this will likely require several Higgs fields and 
various Yukawa couplings, CKM-type explicit flavor vi- 
olations, and further effects. The surplus bound states 
would then become heavy and disappear from the spec- 
trum. The sign problem would emerge, likewise. If it 
can be controlled by the strength of the breaking, how- 
ever, this could provide new insights, a possibility that 
certainly deserves further study in the future. 

For our simulations, which consumed roughly 2000 
core-years, we employ an extension of the available local 



reneralized from QCD 
1. The introduction of 



HMC algorithm for scalars 126 
simulations, for details see [2ll 
temperature and chemical potential proceeds as in QCD. 
Our aim here is a first exploration of the phase diagram 
as a proof of principle. One obstacle is the presence of 
a bulk transition [2J, |26| . 28| , which turns out to persist 
with dynamical fermions [2l|, [23| • To avoid this we chose 
a lattice with at least Nt = 5 time slices in finite tem- 
perature simulations. Simulations of this theory incur 
considerable computational costs so that we restricted 
our lattices to A^^ = 16 points in spatial directions. Al- 
together, we investigated three different sets of lattice 
parameters: 

(a) At zero density we varied the lattice coupling /3 be- 
tween 0.9 and 1.0 in order to control temperature on our 
5 X 16^ and 6 x 16^ lattices. 

(b) At finite chemical potential /i we used /3 = 0.9 on 
both, 6 X 16"^ and 8 x 16^ lattices, and /3 = 1.0 in zero- 
temperature simulations on a 16^ lattice. The hopping 
parameter in the fermion determinant was fixed at k = 
0.15625 in all these thermodynamic simulations. 

(c) For comparison we also varied k from 0.15385 to 
0.15625 on the symmetric 16^ lattice with ^ = and 
/3 in 0.85 to 1.1. glj. 

More results from smaller lattices to assess systematic 
effects will be reported elsewhere. 

We measured three observables to study the phase di- 
agram. One is the Polyakov loop. Unlike QCD, due to 
the trivial center, it is not an order parameter of the 



FIG. 2: The Polyakov loop and the unrenormalized chiral 
condensate at zero chemical potential (left panel) and the 
renormalized chiral condensate (right panel). 

quenched theory [20] , but it nevertheless reflects the cor- 
responding first-order phase transition very well [2J, ^^ . 
In fact, it remains so small in the low-temperature phase 
that it is only possible to determine upper bounds. We 
found this to be true also with dynamical quarks. The 
second observable is the chiral condensate. The quenched 
G2 theory has only one first-order transition at finite tem- 
perature which manifests itself also in the chiral conden- 
sate [2^ , a feature that it shares with QCD and two-color 
QCD. This is in contrast to QCD with adjoint quarks, 
where there is no sign problem either |13| . but where 
separate chiral and deconfinement transitions occur at 
largely different temperatures |30l - [33 | . We normalize the 
chiral condensate to its (/3-dependent) vacuum value, to 
avoid explicit renormalization. The third observable is 
the baryon number density, the derivative of the par- 
tition function with respect to the chemical potential. 
At large chemical potentials the density saturates to a 
temperature- independent value. This is observed also in 
two-color QCD [33. To eliminate the scale, we therefore 
normalize the baryon density to its saturation value. 

Saturation on finite lattices can occur due to the Pauli 
principle. In two-color QCD, this happens when all sites 
are occupied by bosonic baryons such that the lattice has 
the maximal filling of ANf quarks per site as allowed by 
the Pauli principle [l^- Here, we also find the analogous 
saturation at a maximum of 14 quarks per lattice site 
for the two spin states of a single quark fiavor in the 7 
dimenisonal fundamental color representation of the G2 
gauge group. However, the fermion dynamics does not 
fully freeze out because the plaquette average at satura- 
tion density is 0.5647(1), which is smaller than the cor- 
responding value of 0.5724(5) for the pure gauge theory. 

In order to assess the phase diagram it is necessary to 
fix at least a relative scale. Using bound state masses for 
this purpose is rather challenging. We have determined 
the masses of the lowest states in the scalar meson and 
diquark channels on the 16^ lattice, see Figured] These 
include only the connected contributions, however. The 
true scalar meson in our one-fiavor case contains discon- 
nected contributions as well which will cause the domi- 
nant splitting of its mass from that of the light diquark. 
The results show strong systematic effects [21|. More- 
over, the diquarks are most sensitive to the quark mass 
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FIG. 3: Bottom row; raw Polyakov loop data (left), chiral condensate (middle), and normalized baryon density (right) at finite 
quark chemical potential fj, and temperatures T = 0, 98 MeV and 131 MeV. Top row: corresponding numerical derivatives. 



due to their would-bc-Goldstone nature [34|. In order 
to fix the scale from these masses would require a much 
more careful and expensive mapping of the lines of con- 
stant physics in j3 and k. Here we simply use the first 
excited state in the diquark channel as the signal with 
the least sensitivity to the lattice parameters. It varies 
for /3 = 0.9... 1.0 at K = 0.15625 from 2.3t| to I.65I5 in 
lattice units. Requiring that the critical temperature for 
/i = is at 160 MeV then leads to lattice spacings be- 
tween 0.25^4 fm to O.I8II20 fi^^ in tt^c same range for /3, 
with an acceptable systematic uncertainty for this first 
exploratory investigation. 

The results at zero chemical potential are shown in 
Figure [D The normalized chiral condensate only shows a 
rather weak response to the transition observed in the 
Polyakov loop. In fact, after a slight drop the chiral 
condensate starts rising again as temperature is further 
increased. That this behavior is likely to be a lattice 
artifact is seen in particular by comparing the data from 
the 5 X 16"^ and 6 x 16'^ lattices as shown for the same 
range of couplings /3 = 0.9. ..1.0 at k = 0.15625 in Fig. [51 

On one hand, this apparently unphysical behavior of 
the condensate might be due to finite volume effects, be- 
cause the three-dimensional volume shrinks with increas- 
ing temperature when controlled by the lattice coupling 
as it is done here. On the other hand, in the zero temper- 
ature simulations we observe that the meson and diquark 
masses at fixed k = 0.15625 start increasing with /3 again 
just above /3 « 1, c.f. Fig. [T] This might indicate that 
we encounter a transition into an unphysical phase such 
as the Aoki phase [21|. Although we have not carefully 
mapped out the lines of constant physics, for the highest 
temperatures in Fig. [21 which correspond to /3 = 1, we 
might be getting too close to or even into this unphysical 
phase. Both, finite- volume effects and Aoki phase would 



lead to an unphysical behavior as observed, a weakening 
of the chiral transition together with an increase of the 
chiral condensate towards higher temperatures. 

While this implies that even larger and finer lattices 
will be required to reliably quantify the systematic uncer- 
tainties in the high temperature regime, we avoid these 
lattice artifacts in our finite density simulations below as 
much as possible with controlling temperature by vary- 
ing the number of time slices at reasonably small /? val- 
ues. From our present analysis we are confident that the 
lattice parameters (b) used at finite chemical potential, 
as described above, should be reasonably safe. A more 
detailed quantitative analysis, in particular at tempera- 
tures near the ^ = transition in Fig.[2land above, seems 
prohibitively expensive at the moment. 

As a final comment on the ^ = transition before 
we discuss our finite density results we note that, strictly 
speaking, there is no chiral symmetry in the usual sense in 
our one-flavor theory. Rather, the condensate here mea- 
sures the breaking of the extended SU(2) Pauli-Giirsey 
symmetry down to the usual baryon number U(l). Thus, 
our intuition about the usual chiral restoration at the de- 
confincmcnt transition might not be a good guide in this 
case cither. It is quite possible, for example, that the 
2-flavor theory behaves differently because it includes a 
standard chiral SU(2)xSU(2) component in its extended 
flavor SU(4) which is expected to be broken at low tem- 
peratures and which might therefore show a stronger re- 
sponse to the fj. = transition in G2 as well. 

The behavior of Polyakov loop, chiral condensate and 
quark density at finite chemical potential with the scale 
set as explained above and in lattice units are shown 
in Figures [31 and [H At zero temperature the Polyakov 
loop peaks at around 1.5 GeV. Chiral condensate and 
baryon density both show a rather broad transition near 
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FIG. 4: The data from the bottom row of Figure [3] in lattice units. 



the maximuni in the derivative of Polyakov loop. At 
large chemical potential the density saturates, as dis- 
cussed above. This effectively quenches the theory which 
explains the decrease of the Polyakov loop as in two-color 
QCD [i^, [331- It also implies, however, that the regime 
beyond the peak in the Polyakov loop is affected by lat- 
tice artifacts thus limiting the chemical potential to below 
1.5 GeV at T = and to below about 1 GeV at the two 
finite temperatures considered here. This is emphasized 
by the result in lattice units where the saturation is seen 
to be nearly independent of the lattice coupling. 

In the zero temperature results for the quark density 
one should see the onset of Bose-Einstein condensation 
of diquarks when the quark chemical potential reaches 
half the diquark mass. With a mass of the lightest di- 
quark around 300 — 400 MeV here, this so-called silver- 
blaze point [3al is thus expected around fi ~ 150 — 200 
MeV. On smaller 16 x 8'^ lattices, with larger masses and 
with considarably better statistics, one indeed observes 
this onset of diquark condensation in the baryon density 
which increases from zero to a small finite plateau in a 
narrow region [21|. The critical chemical potentials for 
this onset furthermore agree very well with half the di- 
quark masses as extracted from their correlators for a 
wide range of lattice parameters [2l|, [S^ . It is a small 
effect, however. So small that it cannot be seen on the 
scales of the Figures presented here. A more detailed 
study with finer resolution on the present lattices requires 
much higher statistics. Having said that, the onset of 
the large increase in density observed here, for T = at 
/U ~ 1 GeV, is by the same argument unlikely to be due 
to these light pseudo-Goldstone diquarks. The first ex- 
cited diquark states of Fig. [1] should come in at fi around 
400 — 550 MeV and thus still tend to be too small to 
explain this threshold. The intriguing alternative would 
be that we do sec the effect of fermionic baryons here. If 
the quark-gluon hybrids and other baryons are too heavy 
as expected, one would conclude that the density in the 
region around the increase at /i « 1 GeV is dominated 
by three-quark states, the C?2 nucleons. We will find out 
with further baryon spectroscopy in the future. 

Finally, when both temperature and chemical poten- 
tial are non-zero, a gradual shift towards smaller chemi- 
cal potentials along with a broadening of the transition 
is observed. The interpolated phase diagram in Figure [5] 
exhibits a rather polygonal shape but is potentially af- 



fected by systematic effects in the scale setting. Results 
on smaller lattices show the same general features albeit 
with larger systematic errors [21[ . 

In total, the phase diagram of G2 QCD shows a number 
of interesting features which deserve further systematic 
investigation to see whether they are genuinely physi- 
cal, or whether some of them are lattice artifacts. Given 
sufficient computational resources its complete phase di- 
agram can be mapped with any desired precision. In 
particular, the question whether one can have a chiral 
first-order transition rather than a smooth crossover at 
finite density is very significant as it would imply the 
presence of a critical endpoint in a theory with proper- 
ties so close to the ones of QCD, where this is of prime 
importance for future heavy-ion experiments. 

To summarize, we have determined for the first time, 
but still at an exploratory level, the phase diagram of 
G2 QCD, a non-Abelian gauge theory with fermionic 
baryons. This opens new avenues for high density stud- 
ies of the strong force from first-principles calculations. 
In particular, it will be possible to assess the role of the 
Fermi statistics at finite baryon density which is of great 
importance, e.g. for compact stellar objects. One might 
ask at which relative densities quark or hadron equations 
of state are more favorable, and thus whether only neu- 
tron stars or also quark stars can exist. Given the inflow 
of new astrophysical observational data on compact stel- 
lar objects, including major new discoveries, any insight 
on a fundamental level would be timely and important. 
Of course, G2 QCD is not QCD, and such insight can 
only be qualitative. 

It also has to be assessed carefully to what extent the 
particle spectra of the two differ and what effects this has. 
However, these ab-initio results provide a novel approach 
to test investigations of the true QCD phase diagram by 
effective models and functional methods. If G2 QCD is 
continuously deformable to QCD [20, lU, [Sg , this may 
proceed along the lines developed for Yang-Mills theory 
[5|-|7[: Systematically check approximations and model 
assumptions in functional methods along this deforma- 
tion from lattice calculations as long as possible, so that 
the latter can then be used reliably in the final step to- 
wards QCD when the sign problem strikes back. 
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